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Abstract 

We construct a model where cosmological perturbations are analytically solved based on dilatonic 
brane worlds. A bulk scalar field has an exponential potential in the bulk and an exponential 
coupling to the brane tension. The bulk scalar field yields a power-law inflation on the brane. 
The exact background metric can be found including the back-reaction of the scalar field. Then 
exact solutions for cosmological perturbations which properly satisfy the junction conditions on the 
brane are derived. These solutions provide us an interesting model to understand the connection 
between the behavior of cosmological perturbations on the brane and the geometry of the bulk. 
Using these solutions, the behavior of an anisotropic stress induced on the inflationary brane by 
bulk gravitational fields is investigated. 
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I. INTRODUCTION 



The possibility that our universe is a hypersurface (brane) embedded in a higher- 
dimensional spacetime (bulk) has recently attracted much attention [1]. Particularly, a 
model proposed by Randall and Sundrum has attractive features for gravity and cosmology 
[2]- [4]. In their model, a single positive tension brane is embedded in 5-dimensional Anti-de 
Sitter (AdS) spacetime. Ordinary matter fields arc assumed to be confined to the brane and 
only the gravity can propagate in the bulk. A quite interesting feature of their model is that 
a 4-dimensional gravity is recovered in the low energy limit even though the size of the bulk 
is infinite. We need no longer a compactification of the extra dimension. Then cosmological 
consequences of their model have been intensively investigated. In string models, the grav- 
ity enjoys the company of scalar field such as dilaton. Then the extension of Randall and 
Sundrum model to dilatonic brane worlds also attracts much interest [5]-[13]. The dilatonic 
brane world provides us a new models for infiationary brane world which is often called bulk 
infiaton model [14]- [23]. The infiation on the brane is caused by the infiaton in the bulk and 
the bulk spacetime itself is not infiating. It has been shown that this model indeed mimics 
the 4-dimensional infiation model at low energies where the Hubble horizon on the brane 
is sufficiently longer than the curvature radius / in AdS bulk . 

One of the most important quantities which should be clarified in brane world models is 
the cosmological perturbation because it provides us a possibility to test the brane world 
model observationally [24]- [29]. Unfortunately, it is quite difficult to find the solutions for 
cosmological perturbations because we should consistently take into account the perturba- 
tions in the bulk. Recently, Minamitsuji, Himemoto and Sasaki investigated a behavior of 
cosmological perturbations in a model with bulk scalar field in AdS spacetime [22]. They 
used a covariant curvature formalism and showed that 4-dimensional results arc recovered 
at low energies. Thus, to predict signatures specific to the brane world model, we should 
investigate the higher energy effects. However, it is quite difficult to perform the calculations 
at high energies because we should treat completely a 5-dimensional problem. Technically, 
it is necessary to solve complicated coupled partial differential equations to find the behavior 
of perturbations. 

Hence it is eagerly desired to construct a model where we can solve the equations for 
perturbations analytically. Recently, we have developed such a model based on dilatonic 
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brane worlds [20]. The action for this model is given by 

S^l d'xV^, (^i? - ^d^cl>d^cl> - A(0)) - J d'x^A{(l>), (1) 

where is five- dimensional gravitational constant. The potential of the scalar field in the 
bulk and on the brane are taken to be exponential: 

= {j + s) >^le-'^"'^t (2) 
k'A(0) - V2Aoe-^'"^^. (3) 

Here Aq is the energy scale of the potential, b is the dilaton coupling and we defined 

A = 46^-^. (4) 

We assume the Z2 symmetry across the brane. The bulk scalar field (f) acts as an inflaton. For 
5^0, the brane undergoes a power-low inflation. The metric for 5-dimensional spacetime is 
written by separable functions of time and extra-coordinate y. This is technically essential 
because the perturbations can have separable solutions which enables us to derive solutions 
analytically. It is not sufficient to derive general solutions for perturbations in the bulk. We 
should find particular solutions which satisfy the correct boundary conditions at the brane. 
It is in general difficult to find such particular solutions for an expanding brane. But, in 
this model, we can successfully find such solutions analytically. The first aim of this paper 
is to provide exact and analj^ic solutions for cosmological perturbations in this model. 

Using the solutions for perturbations, we can address the primordial fluctuations gener- 
ated during inflation. In a previous paper, we calculated a spectrum of curvature pertur- 
bation TZc by quantizing a canonical variable for the second order action [20]. We found 
that ,even at high energies, the effects of Kalzua-Klein (KK) modes are neghgible at long 
wave-length even though the amplitude of the fluctuation is amplifled. However, in the 
brane worlds, the curvature perturbation alone does not determine the cosmic microwave 
background (CMB) anisotropics. The anisotropic stress induced by bulk gravitational fields 
affects the CMB anisotropics. Indeed, it is this anisotropic stress which gives distinct fea- 
tures in CMB anisotropics in brane worlds from 4-dimensional models [30] [31]. Thus we 
should specify the initial condition for the anisotropic stress generated during inflation as 
well as the curvature perturbations. The second aim of this paper is to investigate the 
behavior of anisotropic stress in this model. 
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The structure of the paper is as follows. In section II, the background spacetime is 
described. In section III, the solutions for cosmological perturbations are derived. Section 
IV is devoted to the investigation of the anisotropic stress on the brane which is generated 
during inflation. In section V, wc summarize the results. In appendix A, 5-dimensional 
Einstein equations for scalar perturbations are shown. In appendix B, the procedure to 
derive the solutions for scalar perturbations is presented. 

II. BACKGROUND 

We first review a background solution [20]. For S — 0, the static brane solution was 
found [8]. The existence of the static brane requires tuning between bulk potential and 
brane tension known as Randall-Sunrum tuning. It has been shown that for A < —2, we 
can avoid the presence of the naked singularity in the bulk and also ensure the trapping 
of the gravity. The reality of the dilaton coupling requires —8/3 < A. For A = 8/3, we 
recover Randall-Sundrum solution. The value of 5, which is not necessarily small, represents 
a deviation from the Randall-Sundrum tuning. This deviation yields an inflation on the 
brane. 

The solution for background spacetime is found as 



cP{t,y) = cP{t) + E{y). 



(5) 



The evolution equation for background metric a and are given by 




-2V2bK(t> 



(6) 



The solution for a{t) and are obtained as 




(7) 
(8) 




Hot = (-i?7?) 3(^+2) . 



3A+8 



where 



3(A + 2) 



3A + 8 
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and 77 is a conformal time defined by 

3A + 8 ^ 3(A+2) 

' 3(A + 2) ° ^ ^ 

We should notice that power-law infiation occurs on the brane for —8/3 < A < —2. Thus 

in the rest of the paper we shall assume —8/3 < A < —2. The solutions for W{y) and S(y) 

can be written as 

e^(f) =7i(y)3(^, e^^^y^ ^n{y)^, (11) 

where 

n{y)^^-l-^smhHy. (12) 

Here we assumed ^ + 5 < 0. At the location of the brane y — yo the solutions should satisfy 
junction conditions; 

dyW{y)\y=y, = _e^(^o)-V26.H(,o)^^^^ 9,S(y)|,=,, = -e^(^°)-^^'^=(^°)6«-^Ao. (13) 
Then the location of the brane is determined by 

sinh//yo=(-l-^j ■ (14) 

It is quite usefuU to note that the above 5-dimensional solution can be obtained by a 
coordinate transformation from the metric 

ds^ = e^Q^'\dz^ - dr^ + Sijdx'dx^), e'^'^(^) = e'V2bQ{z)^ (^5) 

where 

gQW ^ (sinhif|/o)"^(^(if^)3(^^ (16) 



by 



z = —rj smh{Hy), 

T = —r]cosh.{Hy). (17) 

Because the metric Eq. (15) is simple, it is convenient to solve the perturbations in the bulk 
not directly in Eq. (5) but in Eq. (15) for scalar perturbations. 

The background equations on the brane Eqs.(6) can be described by the 4-dimensional 
Brans-Dicke theory with the action 



2Ki 

where 



.V2bK<P 1 ^2t. _ _x2rA +4 1 



<psD^e^''^^: ^5^= ^,lVe (^bd)^-K6— . (19) 

2b^ A ifiBD 



III. COSMOLOGICAL PERTURBATIONS 



Now let us consider cosmological perturbations in this background spacetime. Taking 
appropriate gauge fixing conditions, the perturbed metric and scalar field is give by 

ds"" = e^^^^) [e^^'"^^^*)^?/^ _ de + e2"(*) (% + hij)dx'dx^\ , (20) 

for tensor perturbations, 

ds^ = e^^^y^ (e''^'""^^'^dy^ - dt" + e'"^ (^T^dydx' + 2Sidtdx' + Sijdx'dx^)) , (21) 

for vector perturbations and 

ds^ = e^^^y^ \^e^V2b^Ht)^l + 2N)dy^ + 2A(it(i|/ - (1 + 2$)^^^ + e^"^'\l - 2^)6ijdx'dx^)\ , 
<f> = 0(i) + S(y) + «-^(50, (22) 

for scalar perturbations, where perturbations are decomposed according to the tensorial type 
of perturbations with respect to 3-space metric Sij. Here Si and Tj are transverse vector 
(V*S'i — and V^Tj = 0) and hij{y, t, x) is a transverse and traceless tensor (/i- — 0, V*/iij = 
0) where Vj is the derivative operator on 3-space metric 5^ . Because each type of variables 
obeys an independent closed set of equations in the 5-dimensional Einstein equations, we 
derive the solutions for tensor, vector and scalar perturbations separately. 



A. Tensor perturbations 



The evolution equation for tensor perturbations is simple. The equation for tensor per- 
turbation hij{y,t,x) — h{y,t)e^P^eij is given by 

„2%/26k. 



/i+ (3d + \/26fi:0)/i + e-2V/i =h" + 3W'h', 



(23) 



where eij is a polarization tensor and dot denotes the derivative with respect to t and prime 
denotes the derivative with respect to y. The junction condition for h{y, t) is imposed as 



9yh\y=yQ — 0. 



(24) 



We can use the separation of variables to solve this equation. The solution for hij is given 

by 

hij = / dmd^ph{m,p) fm{y)9m{vy^''eij, (25) 
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where 



fo{y) = 1, 

^0(77) = i-Hr))-^^ (H^y_^{-pri) + c{p, 0)H^'l^{-pr])) , (26) 

V A+2 A+2 / 

A / P'lt]. (coshHyo) \ 

Uy) = {sinhHy)W^ ~ Q^^ (coshi/yo) ^-^+-^'"'^^^^ ' 

Qr^irj) = (-//77)-^ {H^t\-PV) + c{p, m)Hi^\-prj)) , (27) 

where Pg and are associated Legendre functions, H^^ and H^^ are Hunkel functions 
and 

A 



^~ 2(A + 2)' " (A + 2)2- ^^^^ 

The coefficients h{p, m) and c(p, m) are so far arbitrary and these are determined if one 
specifies the initial conditions and boundary conditions in the bulk. The canonical variables 
for the second-order perturbed action for tensor perturbation is given by 

^ = Yk^- (29) 
Then the second-order perturbed five-dimensional action for the tensor perturbation is given 

by 

It should be noted that the modes with < m < —H/ (A -|- 2) is not normahzable. Thus 
the normalizable modes have a mass gap and the continuous massive modes start from 
m = -///(A + 2); 

H , , 

~ A+2 ^ ' 

B. Vector perturbations 

The calculations of the vector perturbations in this model are similar with those given 
by Ref.[32]. In order to solve the Einstein equation for vector perturbations, it is convenient 
to define a variable 

V, = S[-% = Vei, (32) 
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where Cj is the polarization vector. Then the Einstein equations of (0, i) and {i, y) compo- 
nents are given by 

e'^°'{V + (5a - V2bK<j))V) = p^T, (33) 
where we expand the variables by Cje'^^. The component is given by 

S + {3a + V2bK^)S = e-^^'"^'^(r' + 3W'T'). (34) 
We can easily find a master variable from which all perturbations are constructed; 

^ _ ^-3W-3a-V2bK<l)Q/ 
rp _ ^-3W-3a+V2bK<j>Q^ 

If Q satisfies the evolution equation, 

n - (3d - V2bK^)fl + e-^Vl = ^" - ^Wn', (36) 



the 5-dimensional Einstein equation is automatically satisfied. The junction condition for 
perturbations are imposed as 

Thus the master variable should satisfy 

Q|,=,„ = 0. (38) 
The solution that satisfies the junction condition is obtained as 

n{t,x,y)^ jdmd^pn{m,p) Um{y)vm{t)e'P'' (39) 



where 



Voir]) = i-Hrj)^^ { H^l''_^{-pr]) + c(p, Q)H^^^_^{-p7])\ , (40) 

\ A+2 A+2 / 



uM = {sinhHy)^ [ P^^.^Jcosh Hy) - ^^^_^ Q\ {cosh Hy) 



^^+,,(coshi^|/o) 
Q\^JcoshHyo)''-t 

Vmiv) = i-Hri)^^ {Hil\-pr])+c{p,m)H^^\-prj)) , (41) 



where 



A + 4 



(42) 



2(A + 2)' \\H^ (A + 2)2- 

The second order perturbed action is also written by Q as 

= 1/ y ^y^^^3^^-3Wiy)^V2bn<p(t)^-3a{t)^^-2V2bn4>{t)^'2 _ ^2 _ ^-2a{t)p2^2y (43) 

From this action, we can determine the normahzation of the perturbations. It should be 
noted that the 0-mode solution is not normalizable for vector perturbations. 



C. Scalar perturbation 

The scalar perturbations are more complicated than tensor and vector perturbations 
due to the existence of the scalar field in the bulk. Unlike a maximally symmetric bulk 
spacetime, we cannot find a master variable for scalar perturbations and this causes the 
difficulty to solve the perturbations. Fortunately, in our background spacetime, there is a 
simple coordinate system in the bulk, that is Eq. (15). Thus we can use the metric Eq. (15) 
to find general solutions for perturbations in the bulk. In Ref [9] , it was shown that there are 
variables which make the equations for N,A,^,^, and S(f) in the bulk to be diagonahzed. 
Then by performing a coordinate transformation, it is easy to find general solutions for 
perturbations in our background spacetime. However, we should proceed further to derive 
solutions for perturbations. On the branc, the perturbations should satisfy the boundary 
conditions. In terms of the metric perturbations, the boundary conditions are given by 

^'\y=yo = W\N-V2b6(l>)\y=y,, 

d<f>'\y=y, = 3V2bW'{N-V2bS<f>)\y=y„ 

A\y=y, = 0. (44) 

A problem is that if we rewrite these conditions in terms of variables that make the bulk 
equations to be diagonahzed, the boundary conditions become complicated. Indeed, we will 
find that the boundary conditions are not diagonahzed and also they effectively contain 
time derivatives of the variables. Thus unlike vector and tensor perturbations, imposing 
the boundary conditions is not so easy. This is in contrast to the case for the static brane 
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which sits at constant value of the z in the coordinate Eq.(15). For the static brane, the 
boundary conditions for variables which make the equations in the bulk to be diagonalized 
also make the boundary conditions to be diagonalized and also they do not contain the time 
derivative of the variables. Thus, this complexity of the boundary conditions reflects the 
fact that our brane is moving in the coordinate Eq.(15). This movement of the brane causes 
the cosmological expansion on the brane. Hence, it is an essential part of the calculations 
to find particular solutions which satisfy the boundary conditions at the brane. In fact, it 
has been recognized that this is the central part of the calculations of scalar cosmological 
perturbations. In general, the bulk perturbations are not written by separable functions 
with respect to a brane coordinate and a bulk coordinate. Thus there is almost no hope 
to find particular solutions analytically and it has prevented us from understanding the 
behavior of scalar perturbations on the cosmological brane. In our background spacetime, 
the bulk perturbations are obtained analytically by separable functions and it enables us to 
derive the solutions which properly satisfy the boundary conditions on the brane. 

Now, we will describe our procedure to derive solutions. We first solve the perturbation 
in a static coordinate; 



It is possible to find variables which make the equations diagonal [9] (sec Appendix A-2) 

o;^ = r - 2^', 



The evolution equations in the bulk are obtained from 5-dimensional Einstein equation as 



^ e^^(^)((l + 2V)dz^ - (1 + 20)dT^ + IGdzdT + (1 - 2xlj)5ijdx'dx^), 

= + K'^Scj). 



(45) 



ouisf = r — V2bS(f), 



uja = G. 



(46) 



□50;^ = 0, 



□ 50; AT 



2(A + 4) 1 

A + 2 ^ 
2 1 



(47) 
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where 

By performing a coordinate transformation, the evolution equation in our background space- 
time can be derived 

n^Uc = 0, 

□50;^ = 0, 

2(A + 4) 
- A + 2 i-Hrjysu.h'Hy''^' 
2 

= A + 2(-//,)^sinh^//y"^' ^^'^ 



where 

°^-'^J (a^ + AT^^^^^^^^^J-ia^f + AT2;^a;^+^ 1- 
The solutions for a;^ are given by 

a;e = y dmXe(me,p)(sinhi/7/)ra (^P^^.^^^^^(coshi^7/) + Ceg%^.^(^^)(cosh(i/7/)) 
a;^ = 1 c^m^AA^(mv,,p)(sinhi/|/)M (^P^^^.^^^^^(coshi7|/) + Cv,g%^ 

j dmAAfA{mA,p){smhHy)W (^P'^^^,.^^^^^^^.^{cos\i 
X ( ~ - AT2 Hi^(^^) i-pr]), 
UN ^ J dmNAfN{mN,p){sinhHy)^ (p^t^.^(^^^(cosh i/y) + C'ArQ^:f+.^(^^)(cosh(i/y)) 

(51) 



where 



A . . 



' (52) 



2(A + 2)' ^ ^ \^H^ (A + 2)2' 
The perturbations ^, 5(f), N, A and $ are related to -0, (^0, T, G and by a coordinate trans- 
formation; 

^ = ^, (50 = (50, $ = ^' - iV, 

= r cosh^ Hy — sinh^ Hy + G sinh Hy cosh 
A = -{-H-q)^^ (2(r - 0) cosh //y sinh //y + ^(sinh^ i/y + cosh^ Hy)) . (53) 
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Prom Eq.(46), 50, F, G and are written in terms of uji as 

2 

^ ^ ~ 3(A + 4) V2huc), 
4 

(a; AT + v2hu}c + 36^a;^), 



3(A + 4) 

= ^ - r. (54) 

Then the general solutions for the metric perturbations in the bulk are obtained as 
2 

* " ~ 3(A + 4) ~ ~ V2&c^c), 

AT = ^(^(2 + 3sinh^ + + sinh^ //y) 

+ _|_ 4^ (2 + 3 sinh^ -H"y)ci;c + sinh Hy cosh HycuA, 



3A+8 

A = -(-if77)3(A+2) 



• 1 rr 1 rr 2 2(A + 3) 2^26 ^ 

2 smh Hy cosh //y I ^^^iv + ^ ^ ^ + ^TI^''' 



+ (l + 2sinh2i7|/)a;Aj • (55) 

We should impose boundary conditions on the brane. In terms of cui, the junction condi- 
tions become 

oj'c = 0, (56) 
2cosh//yosinhi/yo / 2 2(A + 3) 2^26 \ 

^ ~ l + 2sinh-i/,o VaTI"- + ^T^"^ + ATI'^J 



a;^ - a;^ = ^ _^ ^ -^ coth iJyo [^^^Af + - tanh HyoujAj , (58) 

(59) 



sinh^ Hyooj'j^ + ^1 + '^^^^^^ ^inh^ -H"?/o j i^V' + ^^"^ "^^o Hy^u;'^ = 0, 



A + 4 

where above equations should be evaluated on the brane y — y^. The variables should 
also satisfy the "constraint equations" which do not include the second derivatives of the 
variables with respect to t and y in 5-dimensional Einstein equations, that is, (t, i), {y, i) and 

{y, t) components of Einstein equaitons. Among them, the equation obtained by combining 
(0, i) and {y, i) components of the 5-dimensional Einstein equation (see Appendix A-2 for 
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derivation) 



(1 + sinh^i7y)(x;^ + 



H coth HyujA + 2 sinh Hy cosh Hyu'^ 



{—Hrf) {cosh Hy sinh Hydj^UA + 2(1 + sinh^ Hy)drfOJ^ , (60) 



will be useful to find solutions. Because oja is a variable which is associated with Z2 odd 
variable A, o;^ in the junction conditions would be eliminated using constraint equations. 
Indeed, by projecting Eq.(60) on the brane, we find that a;^(t/o) can be rewritten in terms 
of uja, ^^'ip and uj^p. An important point is that this equation contains the time derivative of 
the variables. Thus effectively, the boundary conditions contain the time derivative of the 
metric perturbations. It should be noted the junction condition for cUc is decoupled from 
other variables. This variable cUc is the canonical variable for the second order action and it 
is directly related to the curvature perturbations on the brane. Thus we do not need to know 
the full solutions for perturbations in deriving the solutions for curvature perturbation. On 
the other hand, in order to derive the anisotropic stress induced by bulk perturbations, we 
should know the solutions for all variables, as we will observe later. This indicates that the 
anisotropic stress on the brane is complicatedly coupled to bulk perturbations compared 
with the curvature perturbation. 

We describe the procedure to derive the solutions which satisfy the boundary conditions 
and constrained equations in Appendix B. In this section we only show the results. The 
solutions are written as the summation of KK modes with mass m; 




(61) 



where rh — m/H . For 0-mode with m = 0, we get 



c^c(O) 




uja{0) 



A + 4 



sinh Hy cosh Hy{—Hrj) ^+'^H 2A+5 (—prj). 



(62) 
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For massive modes with m > —H/ (A + 2) 
V2{A + 2) 



46 



iv - l){smh. Hy)^(^+^) B'^^^.^{cosh. Hy){-Hr)) Hi^{-prj), 



{sinh Hy)^(^{-Hr})-^ --B'^.^Jcosh Hy)Hi,{-pr}) 



1 I iv — 



2A+3 ■ 
A+2 



IV — 



A+1 

A+2 \ n/i 



2\w- 



1 



A+2 

A 



A+3 \ B\^,,X^os]iHy)Hi^^2{-pn) 



A+2 . 



ujN{m) = {sinh Hy)^(^(-Hr]) ^+2 



A + 2 



BniJcosh Hy)Hi,^2(-pv) 



%v 



A+2 . 



^!^t+,.(cosh//y)//,.(-p77) 



2 \iv 



1 



A+2 - 

A 



IV 



A+3 
A+2 . 



ijjA.{m) = (sinhi7y)2(A+2j(_ff^) a+2 
/ 1 



B^^}Aco^\^Ey)E,,{-W]) 



IV 



A+2 . 



IV — 



A+2 . 



IV 



A+2 / 



where 



B^{cosh Hy) = P^(cosh Hy) 



P^tljcoshHyo) 



Q'^VicoshHyo) 



g^(cosh//y), 



(63) 



(64) 



and Ha is the arbitrary combination of Hunkel functions H^^ and Then the solutions 
for metric perturbations are derived as 

J d^'pdmX{m,p)^{m){y,t)e'^'', 
J SpdmN{rh, p)5(t){m) {y, t)e'P'', 



^{y,t,x 
5(f){y,t,x 
N{y,t,x 
A(y,t,x 



d^pdmM{m, p)N{m) (y, t)i 
j d^pdifiM {fh, p)A{m) {y,t)e 
^{y,t,x)-N{y,t,x). 



ipx 



(65) 



where 
^(0) 

5m 



2\/2h [A + 4 i_ 1 i_ 

— — i-Hn) ^+^H 1 i-pri) + i-Hr)) ^+^H 2A+5(-»r?) 

3(A + 4) LA + 3^ " ^" A + 3^ " 

-if 77) ^+■2 H _2A+5{-p'r]) 



3(A + 4) 
7V(0) = \/2650(O), ^(0) = 0, 



-^-t^f-ifr?) ^+2 if 1 (-pr?) - , 
4(A + 3)^ ^'^ A + 3 



(66) 



and 

*(m)(y,t) 



3(A + 4) 



I^-ifry) A+2(sinhii'y)2(A+2y 
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X 



A + 2 



l^--^^B^\lJcoshHy) + (^u - |±|) B\^JcoshHy)^HU-PV) 



.^^-fji/ V^- ATI/ ^2 -1+-^ 2V A + 2A A + 2 



X B\^JcoshHy)\H,,_2{-pn) 



1 



5(l){m){y,t) = l^l-^'?) ^+2(sinh//y)2(A+2) 



A + 2 



X 



^ — S^^tl, fcoshi/y) - 4 6^ - 77^^) B\_^ (cosh Hy)\HU-PV) 



^^-aT2 "^^'"^ \^ 4(A + 2)y -n^'^' ^'j 



A+2/ A+2 



A + 2 y V A + 2 



X B''_,^JcoshHy) \ Hi,_2{-pv)\ , 

"2(A + 2) 



N{fn){y,t) = (--H"r7) ^ (sinh ify) 2(^^+2) 



3(A + 4) 



ii/-^ -5+^^^' \^ 4(A + 2)y -h+i^^ 



+ 



A+2 

1 \ ( 1 



\B^^f^( cosh Hy) 



X (^-f +..(^°«h i/i/) - {^u - (.z. - ,^,,(cosh i/i/)) sinh^ Hy] 



X i7i^_2(-pry)] , 



3A+8 



iv — ^ 



-B^AijcoshHy)H,,{-prj) 



A+2 . 

1 \ r 1 



+ 2sinhHycoshHy \ -— ^ J }^~^^^B^t^^^^{coshHy) 



A+2 

A + 3 /. 2A + 3\ /. A + 1 

IV ■ IV 



A+4V A+2yV A+2 

1 \ /. A + 1 



^^:|+,.(cosh//y)|//,._2(-p77) 



+ (1 -2cosh' Hy) (— TA+s) (^^ - AT2) B'_'-^^iS^oshHy)H,^.,{-m) 
The solutions for massive modes evaluated on the brane become somewhat simple because 



(67) 



the function B^ satisfies 



5!!t^,,(coshi/yo) = 0. (68) 



Using the equations presented in Appendix B, we get 
*(m)(yo,i) = i(-//77)-sT2(sinh//yo)™5!i+,,(cosh//yo) 
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Hiu-2{-PV) 



$(m)(|/o,t) 



{-Hri)~^ (sinh Hyo) i 




(cosh Hyo) 



S(f){m){yo,t) 




(cosh Hyo) 



Hiu{-pr)) + 



V2b f \ 
A + 2 l^i^ 



1 



A+3 



A+2 



) 



Hiu-2{-pr]) ■ 



(69) 



These solutions are first main results of this paper. They provide us the connection between 
the behavior of the perturbations on the brane and the perturbations in the bulk. 

The second order action for scalar perturbations is written in terms of the canonical 
variable cUc, 



This can be verified using the result for the metric Eq.(15) because ujc does not change by 
the coordinate transformation. This action is the same as the second order action for tensor 
perturbations. Then, the massive modes with < m < —H/ (A + 2) are not normalizable. 
Thus there is also mass gap for the scalar perturbations. 

IV. PRIMORDIAL FLUCTUATIONS IN THE BULK INFLATON MODEL 

In the previous section, we obtained the classical solutions for cosmological perturbations. 
These perturbations properly satisfy the boundary conditions at the brane. However, the 
boundary conditions on the brane alone do not fix the solutions completely. There remains 
a freedom to choose the "weight" M{ih,p) in the summation of KK modes Eq.(61). These 
coefficients arc fixed once one more boundary condition in the bulk is specified. Because the 
brane is infiating, it is natural to specify the boundary conditions for the perturbations by 
quantum theory in the same way as the usual 4-dimensional infiationary model. We have 
already derived the second order 5-dimensional action for perturbations, the quantization 
can be done within the full 5-dimensional theory. 



5S^^^ = / dydtd- 





(70) 
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In the previous paper, we have already carried out the quantization of scalar and tensor 
perturbations. It was shown that the KK modes are well suppressed at large scales even 
if the energy scale of the inflation is sufficiently higher than the scale of the bulk. More 

precisely, the bulk curvature scale and the Hubble constant on the brane are determined by 
the bulk potential and the deviation from the RS tuning, respectively. Thus their ratio, 

S 



r — 



(71) 



A/8 + 5 ' 

characterize the behavior of perturbations. If r is large, we expect the five dimensional 

nature of the perturbations become important. However, we showed that due to the mass 
gap in the KK spectrum, the massive KK modes are hardly excited. Thus, at large scales, 
the behavior of tensor perturbation and curvature perturbation defined by 

7^e = -jOUc, (72) 

<p 

are essentially four- dimensional except for the amplitude of the perturbations. 

Hence, we might expect that this model cannot be distinguish from the usual four- 
dimensional inflationary model. However, in the brane world, the curvature perturbation 
TZc alone does not determine CMB anisotropics. The anisotropic stress Sns induced by bulk 
gravitational fields also affects the CMB anisotropics. The anisotropic stress is measured by 
the difference between $ and ^; 

(*-$)!,=,„ = «:V"57r,. (73) 

where — kXq. Then we should also determine the initial condition for Sirg during the 
inflation. Prom the five-dimensional Einstein equation we find that Stte is related to N; 

N\y=yo = i^^e'"'5n£. (74) 

This implies that it is not sufficient to determine the behavior of the canonical variable ujc 
but we need the solutions for all cui. As already observed, the boundary conditions for cui 
except for lOc are complicated and this reflects the fact that the brane is " moving" . Thus 
we expect that the anisotropic stress can have a distinguishable feature which the curvature 
perturbation Tic does not possess. Because we can derive the solutions for N, it is possible 
to investigate the behavior of 67T£. 
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A. Behavior of the canonical variable ujc 



We first review the quantization of canonical variable Uc- The second order action for ujc 
is nothing but the action for a 5-dimensional massless scalar field. Then the quantization is 
easily carried out. The canonical variable uJc can be expressed as 

x,z)^ j dmd^p [aprhOrh{y)Xm{ty^'' + (h.c.)] . (75) 
Here apfh is the annihilation operator and satisfies the following commutation relation, 

The modes functions are given by 

Oo(y) 



6{p — 'p')5{m — fh'). 



(76) 



1 / A\ 2{a"+2) r /-OO^ ,^2 



Orniy) = 



V2 

Fh 



(77) 



and 



A \ 2{A + 2) 



Xo(^) 



(ier + ICr)-^(sinh//y)mS%^.^(cosh//y), (78) 



2 A+2 



y^H-H-Hrjr^^e-'^Hil\-pn\ 



where 



2(A + 2)' 



(A + 2)2' 



' A+l 
^ A+2 



A , 1 . s , . __ . vre . s 



r(f±i-i^) Qn^^.^(coshi/yo) n^ + iy) 

Now the spectrum of the KK modes //{fh, p) is determined as 



(79) 
(80) 

(81) 
(82) 

(83) 



2tt 

Af{0,p) = K^H- 



A \ 2(A+2) 



(sinh Hy') a+2 dy' 



yo 



J\f{m,p) — —K 



A + 2\iu-lJ\ 85 J 



(84) 



The ratio of massive modes and massless mode increases with r. But the ratio becomes 
constant for large r. In the previous paper, it was shown that this is caused by the existence 



18 



of mass gap. And also, the amplitude of massive modes are dumped after the horizon 
crossing. Thus the spectrum of the massive modes is blue tilted, so the contribution from 
massive modes becomes negligible at large scales. We should note that the integration over 
fh logarithmically diverges. Thus we need some regularization scheme. 



B. Behavior of anisotropic stress 



Now we turn to the anisotropic stress 



-yo- 



(85) 



First let us consider the 0-mode. The 0-mode solution satisfies 



3(A + 4) 



-(-^^)-^^i'U' 



A+2 



-pri) 



(86) 



A + 3 

As mentioned in section II, the effective theory for background spacetime is given by the 
BD theory. In the BD theory the correspondent equation is given by 



^BD 



V2bS(p. 



(87) 



As expected, the 0-mode solution can be described by the BD theory including anisotropic 
stress. At late times —pr] — > 0, A^(0) behaves as A^(0) = const. 
The massive modes also contribute to the anisotropic stress; 



X 



dmJ\f{m,p){smhHyo)^(^+^B'^i_^.^{coshHyo){-Hr]) a+2 

3A+5 



A+2 



-pv) 



IV — 



A+3 
A+2 



H, 



(1) 

iv-2 



-pv) 



(88) 



At the horizon crossing —pr] = 1, the ratio of the amplitude of massive mode to 0-mode 
modes has similar feature with the curvature perturbation. The ratio increases with r, but 
it becomes constant for large r due to the mass gap. Note that the term proportional to 
Hil-2i~PV) does not give an additional divergence in the integration over rh. 

However, the subsequent evolution of the anisotropic stress is quite different from the 
curvature perturbation. After the horizon crossing, the 0-mode remains constant. On the 
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other hand, the massive modes increase as (— p?y)-(2^+5)/(^+2) qc e~^(2^+5)a/2 ^-^e term 

proportional to Hi^-2{—pv)- Thus if A < — |, the massive modes will dominate over 0-mode 
and it seems to leave significant consequences in the inflationary brane world. Indeed, the 

massive modes on the brane Eqs.(69) grow for {—pr]) — > 0. Then one might worry that this 
indicates the gravitational instability of the spacetime. However, the physical amplitude of 
the anisotropic stress is given by 

due oc e-^"iV|,=,„ oc e-(6^+i9)a/2^ ^gg^ 

which always decreases with time for — 2 > A > —8/3. The same situation occurs in the 
analysis of the radion in Randall-Sundrum de Sitter two branes. Let us consider two de Sitter 
branes in AdS^ spacetime. By imposing a fine tuning on the tensions of two branes, the 
distance between two branes, the radion, becomes constant. However, if one considers the 
perturbation of the radion, the radion has negative mass squared. In ref [33], the effect of the 
quantum radion was investigated. They found that the metric perturbations in Longitudinal 
gauge grow due to the instability of the radion. However, the physical amplitude of the 
anisotropic stress itself decays. The resolution is that the Longitudinal gauge is not really 
a good gauge. It is possible to find a gauge where all perturbations do not grow. Thus the 
growth of the metric perturbations does not directly imply the instability of the spacetime. 
In our case, the same arguments should be applied. In general, the curvature perturbation 
TZc is the measure of the linear perturbation amplitude. In our background the curvature 
perturbation does not show the instability. Thus the growth of the metric perturbation is 
the artifact of the bad choice of the gauge. 

Let us explicitly show that we can find a suitable gauge where all perturbations are not 
growing. Let us consider a 4-dimensional gauge transformation 

rj^rj-e", x' ^ x' - e'\ (90) 

By choosing = e^, the metric and the scalar field are transformed as 

ds'^ = e^" (-(1 + 2^)d7f + ((1 - 2^)5ij + E^ij) dx'dx^) , 
(f0 = 50 - K(9^0)e'', (91) 



where 



l> = $ - 9^6" - (a^a)e''. 
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^ = * + {dr,a)e'', 

E = -e. (92) 

It is a non-trivial problem to find that simultaneously eliminates the growing part of the 
$(m), ^{m) and S(f){m). In this case, it is possible to find such a solution. By taking 

e\rh) = --L('__L_'j {smhHyo)mB'^,^JcoshHyo){-Hrj)^^ 

^ {H\Z2{-Pn) + H\l\-pn)) , (93) 

the resultant metric and the scalar field become 

= ^(-Tzi) (smhi/yo)^S%^.^(coshi/yo)(-i/^)-^i/SVp^), 
^("^) = .'''"a+s I (sinhi/yo)™S%+,,(coshi/yo)(-i^77)"^/^S^(-p77), 
5<t>{m) = -_(A + 2) Ui/ 



X (sinh ii-yo) +i^(cosh Hyo){-Hr]) H\^{-pq) , 



(94) 



At late times —pr] 0, E behaves as 

E oc (-77)-^ oc e-^"/2 ^ 0. (95) 

Thus we can find a gauge where all perturbations remain small. The anisotropic shear E 
induced by massive modes decreases in an infiationary background A < — 2. Because the 
amphtude of E at the horizon crossing is suppressed due to the mass gap, the effect of the 
massive modes is always negligible. 

In summary, in the anisotropic stress, the contribution of the massive modes will dom- 
inate over 0-mode and this causes the growth of metric perturbations in the Longitudinal 
gauge on the brane. However, we should carefully choose the gauge in evaluating the effect 
of these massive modes on metric perturbations. It is possible to find a "good" gauge where 
all perturbations remain small. In this gauge, the contribution of massive modes are not 
strong enough to cause the instability on the brane and massive modes do not leave signif- 
icant consequences. So we conclude that, in our bulk infiaton model, the long-wavelength 
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perturbations are quite similar with the perturbations in the BD theory described by the 
0-mode solution even for the high energy inflation. 

V. CONCLUSION 

In this paper, we derived the exact and analytic solutions for cosmological perturbations 
in dilatonic brane worlds. We used a background spacetime where the brane undergoes a 
power-law expansion due to the bulk scalar field. The effective theory on the brane is given 
by a Brans-Dicke theory. The interesting feature of this model is that we can derive the 
exact background solutions including the back-reaction of the bulk scalar field. Moreover the 
spacetime metric is separable with respect to the brane coordinate and the bulk coordinate. 
Then it is possible to solve the cosmological perturbations analytically. 

Scalar perturbations are quite complicated because of the existence of the bulk scalar 
field. We can find variables which make the equations in the bulk to be diagonalized. But 
the boundary conditions for these variables are not diagonalized and also they effectively 
contain the time derivatives of the variables. The exception is the canonical variable cUc of 
the action which is related to the curvature perturbation Tic on the brane. The evolution 
equation for uj^ and the junction condition on the brane are decoupled from other variables. 
However, if one wants to derive the solutions for all metric perturbations, we should solve the 
complicated boundary conditions. Because this complexity is caused by the expansion of the 
brane, the derivation of the solutions for this problem is a central part of the calculations 
of cosmological perturbations in the brane world. This difficulty has prevented us from 
understanding the behavior of scalar perturbations on the brane. In this background, it 
is possible to derive the solution analytically. Then we found the solutions for all metric 
perturbations which properly satisfy the junction conditions on the brane. 

As an application, we have investigated the behavior of the anisotropic stress on the 
brane induced by the bulk perturbations. We used the quantum theory for the 5-dimensional 
perturbations to determine the amplitude of the perturbations. As was shown in our previous 
paper, the massive KK modes do not significantly contribute to the curvature perturbations 
even in the high energy inflation where the Hubble scale on the brane is larger than the 
curvature scale in the bulk. In the anisotropic stress, the contribution of the massive modes 
is suppressed at the horizon crossing. Remarkably, however, the subsequent evolution of 
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the anisotropic stress is quite different from the curvature perturbations where massive 
modes rapidly decays. In the anisotropic stress, the contribution of massive modes seems to 
dominate over 0-mode after the horizon crossing for A < —5/2. The difference comes from 
the junction conditions. The junction conditions effectively include the time derivative of 
the variables except for the junction condition for curvature perturbation. This causes the 
growth of the massive modes. However, it also causes the growth of metric perturbations 
in the Longitudinal gauge. Thus a careful choice of the gauge was required to discuss the 
effects of these massive modes. We found a suitable gauge where all perturbations remain 
small. In this gauge, there is an anisotropic shear which describes the contribution of the 
massive modes. It was shown that the anisotropic shear decays in our spacetime. Thus, the 
contribution of massive modes are not strong enough to cause the instability on the brane 
and massive modes do not leave significant consequences. We concluded that, in our bulk 
inflaton model, the perturbations are quite similar with the perturbations in the BD theory 
described by the 0-mode solution at large scales even for the high energy inflation. 

Our analysis indicates that the behavior of anisotropic stress is quite complicated com- 
pared with the curvature perturbation. The behavior of the curvature perturbation can be 
determined by partially solving the perturbations. But in order to determine the anisotropic 
stress, we needed to derive full solutions for perturbations. This is indeed the generic feature 
of the brane world cosmological perturbations. For example, in the Randall-Sundrum model, 
the behavior of the curvature perturbation is determined only by the conservation of the 
energy-momentum tensor on the brane at large scales. But the anisotropic stress can be de- 
termined only if the gravitational field in the bulk is completely specified. So far the analysis 
of this anisotropic stress is very limited due to the difficulty of solving the full 5-dimensional 
perturbations. Our solutions would provide an interesting toy model for the investigation 
about the relation between the behavior the anisotropic stress on the brane and the bulk 
gravitational field. In this paper, we determine the boundary condition for perturbations in 
the bulk by quantum theory. But it is also possible to consider other boundary conditions. 
For example, it might be interesting to consider the boundary condition which allows the 
existence of dark radiation on the brane. Of course, in the inflationary background, the 
dark radiation does not play a role, but it is certainly important to fully understand the 
relation between the geometry of the bulk and the behavior of anisotropic stress, because it 
plays an essential role in the calculation of the CMB anisotropics in the brane worlds. For 
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this purpose, it may be useful to re-derive our results using a covariant curvature formalism 
because the geometrical meanings is manifest in this formulation. We will report these issues 
in the near future [34]. 



APPENDIX A: 5-DIMENSIONAL EINSTEIN EQUATION 
1. 5-dimensional Einstein equation 

For convenience, we present the 5-dimensional Einstein equation for the scalar perturba- 
tions of the metric 

{t, t) component 



- Se-^/^aN + 3e-^^(2d + 7)^ - Se'^^a'N' - 36"^^ (4a' - 7')^' + Se'^'-'dA' 

- 6e-^^ [a" + 2(a'f - a'V] N + 6e-^^a(a + 7)$ + Se'^'^ [a' + a{3a' + (3' - 7')] A 

= «:2(e-2/3($02 _ 050) _ A'50 + e-^''{N(f)'^ - <j)'5<j)')). (A2) 



{t, i) component 



1 



— e 



-2/3/ 



a - 7)Ar - e-2^(2d + 7)$ - -e-^^(a' + 3/5' - 7')^ 



(j)5(j). 



(A3) 



{t, y) component 



3e-2/3^' - 3e-2/^a'iV - 3e-2^(a' - /5')^ - 3e-2^d$' + ^e-^^VM 
6e-2^(Q;' + aa' - a(5' - a'^)^ - 36-'^ \a" + {a'f - a' {(5' + 7') 



A4>'^ + 200'$ - 050' - 0'50 



(A4) 



(i, i) component 
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- e-^'^{2a - /3 + 2^)N + 2e-^^{3a - /3 + 7)* + e-^^{2a + 7)6 + e-^'^{2a' + 2/3' - ^')A 

- e-^^2a' + (3')N' - 2e-'^(3a' + - 7')^' + e-^^2a' + 2(3' - 7')$' + e-2^(2d + 



2e-2^ \2a" + p" + 3(a')' + 2a'/3' - 2aV + {P'f - P'i 



N 



+ 26'"^^ [2ci + 7 + ?,d? - 2aP + 2a^ - /37 + 7^ 



+ e 



-27 



4q;' + 2/3' + 2d;(3a' + /5' + 7') + l3{2a' + 2/3' - 7') - 7(2a' + /?') 



A 



— e 



-2/3 



(i j) component 



{i, y) component 



TV - ^ + $ = 0. 



1 ; 



(A5) 



(A6) 



2 xi,'-^' - A+ {2a' + p')N + {a' - /?')$ - -{a + p + ^)A 

= k20'50. 

(y, y) component 

3 - 2e-2"V2* + e'^^V^* + 3e-2^(4d - /3)^ + 3e-2^d6 + 'ie-'^^olA 

- ?,e-^\2a' + /3')*' + 3e-'^Q;'$' - 6e-'^Q;'(Q;' + I3')N 



(A7) 



+ 6e-2^(a + 2a^ - a/3)$ + 36"=^'^ \a' + a'(3d + /^ - 7) 



K 



- $02) - A'Scj) - e-^^{(f)'^N - (f)'S(f)') 



(A8) 



2. Equations for 

In order to derive the evolution equations for uji in the bulk, it is easy to use the coordinate 
Eq.(15). By combining the Einstein equations and the equation of motion for the scalar field, 
we first get the evolution equations for metric perturbations ip, (p, T, G and 50; 

2 _^dA 



92V' + 3(a,Q)9,V-p'V'-92V' = -2(92Q)r-2«(9,g)(9,0)50 + 



-K 



dcj) 



dlT-r2>{d,Q)dj:-p'T-dlT = -{dlQ-r2>{d,Qf-K\dAf)T 

1 
3 



K {(d,Q){dA) - 9^0) 50 - ^«-'^e2«50. 



dlS(f) + 3{d,Q)dM-p'^S(f)-d^S(f) = 2K{dl(P)r + 2K\d,(f)fS(f) + k 
dlG + 3(a,Q)a,G' - p^G - ^^G = -^dlQ)G 



#2^ 



:2«50, 



(A9) 
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These equations can be diagonalized using Ui defined in Eqs. (46) [9]. Using 

we get the evolution equations for uj^. 

The constraint equations are also easy to be derived using the metric Eq.(15). The (r, i) 
component of Einstein equation is given by 

- \{d,G + ?,{d,Q)G) - 2driJ + drT = 0. (All) 

Rewriting these equations by cu,, we get 

- \{d,u:A + 3(9,Q)a;^) + d-^u:^ = 0. (A12) 

Then performing the coordinate transformation, we get Eq.(60). The remaining two con- 
straint equations can be derived in the same way or can be obtained directly in our spacetime. 



APPENDIX B: DERIVATION OF SOLUTIONS FOR SCALAR PERTURBA- 
TIONS 

1. 0-mode 

In order to derive the solution the formula for the derivative of Hunkel functions is usefuU; 

^ {-Hr])~'^H__^{-pri) = -p(-H7])~^H A+aJ-pr]) 

1 1 1 

-- p{—Hri) ^+2H A+3 (—pr}) 



d_ 
dr) 



[—Ht]) ^+'^H 2A+5 (—prj) 

A+2 



A+2 



2(A + 3) , , A+3 , , 

+ \^H{-Hr))-^^H_._A^{-pr)). 



(Bl) 



Let us find the solution for cui with 

u;,^(-Hnr^^H_^{-pn). (B2) 



A+2 



First we use the constraint equation Eq.(60). Because this equation contains the first deriva- 
tive with respect to time, the solution for cui necessarily includes H 2A+5 (—prj) because it is 



A+2 



only the choice that can eliminate a+3 (— ^77) which arises when we take the time derivative 
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of H i_{—pr]). Then we assume that cui contains H i_{—pr]) and H 2A+5 (—prj). Using 

the fact that cui should satisfy the evolution equation in the bulk, the y-dependence of the 
variables is automatically determined. Then substituting the ansatz for the variables into 

the junction conditions, wc can determine all coefficients except for an over-all normaliza- 
tion. It is a non-trivial check that these solutions indeed satisfy three constraint equations. 
It is easy to verify that these solutions indeed satisfy the constraint equations. 



2. Massive modes 

For massive modes, the following formula is useful; 

d_ 
dr] 



{-Hrj) Hi^(-pr])\ = -p(-Hr]) Hi^_i(-pr]) 

1 



+ H 



d_ 
dr) 



.A + 2 
/2A + 5 



+ w){-Hr]) ^+^Hi^{-pr]). 



+ H 



V A + 2 



A+3 

H-q) ^+^H^2+iu{-pv)- 



(B3) 



and 

d 
dy 

X 

d 
dy 
d_ 
dy 



(sinh Hy) ^^^+^) 5^+^ (cosh Hy) 



-coi\iHy^^B>;+\cos\iHy) + 1^(3 
(sinh 2(^S^(cosh Hy) 



i?(sinh Hy) 2(^^+2) 
A + 1 \ / 



2(A + 2)^ 



/3 



3A 



2(A + 2), 



B^0^\cos\iHy) 



(sinh Hy) (cosh Hy) , 



(sinh //y) 51^5^+1 (cosh Hy) = //(sinh Hy)^(^ [coth HyB^'^\cosh Hy) + S^+^(cosh Hy)] , 



X 



//(sinh //y)2(A+2) 



(B4) 



As for the 0-mode, the ansatz for the solutions are determined so that they can satisfy the 
constraint equation Eq.(60). We asssume 



= i-Hrj)-^^ (sinh Hy)mC,B^^,_^.^Hi,(-prj) 
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X 



C^B\^.^{cos\iHy)H,,{-pi^) + D^B^_, ^Jcosh Hy)Hi,_2 



+11/ ^ 



ua = {-Hrf) A+2(smhify)2(A+2T 



X 



a;jv = {-Hr]) A+2 (sinh i^y) 2('^^+2) 



X 



Civ5^:t+,,(coshi/y)//,,(-p?7) + DMBniJcoshHy)H, 



iu-2 



(B5) 



where B^ is given by Eq.(64). First let us consider the constraint equation Eq.(60). Because 
this equation only contains cu^ and cua, it is easy to determine the coefficients using the 
formula for derivatives. We get 

A+l\ 



If 1 \ 1 / 2A + 3\ Ai/- 

Cjh — — iu — Ca, Dji, — - [iv Da-, Da — — 

^ 2V A + 2y ' ^2V A + 2y ' 



Next we use the junction conditions. From Eq.(58), Dat is determined as 



A+2 ^ 
A+2 . 



D 



1 



1 



AT 



2 \iv — 



A+3 
A+2 . 



From Eq.(57), C^r is given by 



A + 2 
Cn — — — Ca- 



(B6) 



(B7) 



(B8) 



We should note that at this time, the problem becomes non-trivial because the equation 
should be satisfied by the coefficients which have already determined. The point is that, on 



the brane, -B^tl- (coshif|/o) = 0. Using this fact, we can show that 



sinh if I/O cosh HyoB^^'.-icoah Hyo) 



+11/ ^ 



2(iy-l) \" A + 2 
A + 1 



A+ 1\ / 2A + 5 
— I r — 



A + 2 



B% Jcosh Hyo) 



IV 



A + 2 



cosh2//yoS^:t+,,(coshii|/o). 



(B9) 



Using this relation, it is shown that the junction condition Eq.(58) can be satisfied. Finally, 
we use Eq.(57). We get 



C.= -f(A + 2)(i.-^)(i.-l). 



(BIO) 



Here we again used Eq.(B9) in order to show that the junction condition is satisfied. Now all 
coefficients are determined except for an over-all normalization. It remains a task to verify 
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that these solutions satisfy the remaining two constraint equations. In order to show that, 
we use the formula for associate Legendre function and Hankel function 



Hp.iiz) + Hp+i{z) = 2f3z-'Hf,{z). (Bll) 

Then after long calculations, it is possible to show that the above solution indeed satisfy the 
constraint equations. 

In order to derive the solutions for metric perturbations on the brane, we used Eq. (Bll) 
and the following equations; 

^ {'''-^^){^(i^-^)^^^'Hyo+(i^u-^±^^^B^^^^^^^ (B12) 

B%Jcosh Hyo) = (iu - ^) { (iu + ^) +2{zu- 1) sinh^ Hyo] B^^Jcosh Hy,), 

(B13) 

which can be derived using B'^^ (cosh Hyo) — 0. 
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